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Following our previous work [S. Murakami, N. Nagaosa, S. C. Zhang, Science 301, 1348 (2003)] 
on the dissipationless quantum spin current, we present an exact quantum mechanical calculation 
of this novel effect based on the linear response theory and the Kubo formula. We show that it is 
possible to define an exactly conserved spin current, even in the presence of the spin-orbit coupling 
in the Luttinger Hamiltonian of p-type semiconductors. The light- and the heavy-hole bands form 
two Kramers doublets, and an SU(2) non-abelian gauge field acts naturally on each of the doublets. 
This quantum holonomy gives rise to a monopole structure in momentum space, whose curvature 
tensor directly leads to the novel dissipationless spin Hall effect, i.e., a transverse spin current is 
generated by an electric field. The result obtained in the current work gives a quantum correction 
to the spin current obtained in the previous semiclassical approximation. 

PACS numbers: 73.43.-f,72.25.Dc,72.25.Hg,85.75.-d 



I. INTRODUCTION 

Spintronics, the science and technology of manipulat- 
ing the spin of the electron for building integrated in- 
formation processing and storage devices, showed great 
promise^. Spintronics devices also promises to access the 
intrinsic quantum regime of transport, paving the path 
towards quantum computing. However, many challenges 
remain in this exciting quest. Among them, purely elec- 
tric and dissipationless manipulation of the electron spin 
and its quantum transport is one of the most important 
goals of quantum spintronics. 

In our previous work 2 ", we discovered a basic law of 
spintronics, which relates the spin current and the electric 
field by the response equation 

fj = a s eijkE k (1) 

where j* is the current of the i-th component of the spin 
along the direction j and e^-fc is the totally antisymmetric 
tensor in three dimensions. Sinova et al£ found a simi- 
lar effect in the two-dimensional n-type semiconductors 
with Rashba coupling. This law is similar to Ohm's law 
in electronics, and the spin conductivity a s has the di- 
mension of the electric charge e divided by the scale of 
length. However, unlike the Ohm's law, this fundamen- 
tal response equation describes a purely topological and 
dissipationless spin current. It is important to note here 
that the spin current is even under the time-reversal op- 
eration T. When the direction of the arrow of time is re- 
versed, both the direction of the current and the spin are 
reversed and the spin current remains unchanged. Since 
both the spin current and the electric field in Eq. (JTJ are 
even under time reversal T, the transport coefficient a s is 
called "reactive" and can be purely non-dissipative. This 
is in sharp contrast to the Ohm's law 

ji = vEi (2) 



relating the charge current ji to the electric field. In 
this case, the charge current ji changes sign under time- 
reversal T, while the electric field is even under T. Since 
the Ohm's law relates quantities of different symmetries 
under time reversal T, the charge conductivity a breaks 
the time-reveral symmetry and describes the inevitable 
joule heating and dissipation. Quantum Hall current, 
which is transverse to the electric field and dissipation- 
less, has the feature similar to Eq. QJ, but the time- 
reversal symmetry is compensated by the external mag- 
netic field. Dissipationless current without time-reversal 
symmetry breaking is extremely important and funda- 
mental in solid-state physics, the most celebrated exam- 
ple of which is the superconducting current. It is de- 
scribed by the London equation, 

P'^ a tp ldAi r*\ 
Ji = A , E t = — 3 

mc c at 

where the current ji is related to the vector potential Ai 
instead of the electric field. In the London equation, both 
ji and Ai are odd under time reversal T, therefore, the 
transport coefficient p s , also called the superfluid density, 
describes the reversible and dissipationless flow of the 
supercurrent. 

In summary, the dissipationless spin current discovered 
in Ref. shares some basic features with the supercon- 
ducting current and the quantum Hall edge current, in 
the sense that, 1) the spin Hall conductivity a s is a dis- 
sipationless or reactive transport coefficient, even under 
the time-reversal operation T; 2) the spin Hall conduc- 
tivity cr s can be expressed as an integral over all states 
below the fermi energy, not only over states in the vicin- 
ity of the fermi energy as in most dissipative transport 
coefficients. Furthermore, just like the case of the quan- 
tum Hall effect^, the contribution of each state to the 
spin Hall conductivity a s can be expressed entirely in 
terms of the curvature of a gauge field in momentum 
space, which in our case is non-abelian. The dissipation- 
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less spin current is induced by the electric field through 
the spin-orbit coupling, whose characteristic energy scale 
exceeds the room temperature in many semiconducting 
materials. 

Electronic structure of semiconductors with diamond 
structure (e.g. Si, Ge) and zincblende structure (e.g. 
GaAs, InSb) are well understood in terms of the k ■ p 
perturbation theory. The top of the valence bands are at 
k = 0, i.e., T-point. They consist of the three p-orbitals 
p x ,p y ,Pz with spin up and down. In the presence of the 
relativistic spin-orbit coupling, these 6 states are split 
into four- fold degenerate S — 3/2 states and two- fold de- 
generate S — 1/2 state. Here S denotes the total angular 
momentum of the atomic orbital, obtained through the 
coupling of the orbital angular momentum L and the spin 
angular momentum s. The second order perturbation in 
the k • p results in the effective Hamiltonian near k = 0, 
which is called Luttinger Hamiltonian^: 



H =^cj 4ik i? (W ,(k)c I/)k 



#^(k) 
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(7i + ^72)£ 2 - 2 72 (k-S) 



(4) 
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where k = {k x , k v , k z ), S = (S x , S y , S z ), and k = |k|. 
The explicit form of the matrices S l (i = 1, 2, 3) is given 
in Appendix ^ For simplicity, we have put 72 = 73 in 
the original Luttinger Hamiltonian; most of subsequent 
discussions are also applicable to more general cases with 

72 ^ 73- 

On the other hand, the conduction bands are made out 
of the s-orbital and hence doubly degenerate. When we 
neglect a small effect due to broken inversion symmetry, 
this degeneracy is not lifted due to the Kramers theorem. 
Therefore the effect of the spin-orbit interaction is small 
in the conduction band, although the Rashba eSecil^^i^ 
is induced by the electric field near e.g. the interface 
structure. The spin Hall current in the Rashba system 
has recently been discussed by Sinova et alA. These au- 
thors showed that dissipationless and intrinsic spin Hall 
current can take an universal value in this system. The 
position of the conduction band minima depends on the 
material. For example, they are located at general points 
along the axis between the T- and the A-points in Si, 
while they at at the L-points in Ge. We will focus on the 
valence bands below because of the intrinsically strong 
spin-orbit interaction. 

Although the band structure of semiconductors with 
spin-orbit coupling has been understood for many years, 
only recently has it been recognized that the gauge field 
and its curvature in the momentum space made out of the 
Bloch wavefunction play important roles in the transport 
properties of electrons in solids. The gauge field is defined 
in terms of the Bloch state Ink) as 



A n i (k) 



nk 



d 



nk 



(5) 



where n is a band index. It represents the inner 
product of the two Bloch wavefunctions infinitesimally 
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FIG. 1: Schematic picture of the band structure of GaAs 
near k = 0. CB means the conduction band, HH the heavy- 
hole band, LH the light-hole band, and SO the split-off band, 
respectively. When the small inversion symmetry breaking 
is neglected, all of them are doubly degenerate. The Fermi 
energy shown in the figure corresponds to the hole density 
~ 10 19 cm~ 3 . The splitting A = 0.34eV at k = between the 
LH,HH and SO bands are due to the spin-orbit interaction. 
The HH and LH bands are degenerate at k = 0, giving rise 
to a monopole in the gauge field as discussed in Section ITTTl 



separated in k-spacc. This gauge field A n i(k) de- 
scribes topological structure of Bloch wavefunctions in 
the momentum spaced, and plays an important role 
in transport p r0 pertie9ii*i2iiiii2iiiiii£iiLI2ii2i and in mag- 
netic superconductors^. In particular, this gauge field is 
related to the transverse conductivity a xy as 



-E 

h H 

n.k 



n F (e n (k))B nz (k), 



(6) 



where B nz (k) = F n>xy (k) = - is the z- 

component of the field strength made from A n i(k), and 
^F(en(k)) is the Fermi distribution of ff the n-th band 
with energy e n (k). This formula^ is the foundation of 
the integer quantum Hall effect (QHE), and further ap- 
plied to the anomalous Hall effect (AHE) in ferromag- 
netic metala 11 i 12 d 3 i 14 d 5 i 16 d 7 d 8 i 19 . Especially in the mag- 
netic semiconductors (Ga,Mn)As2i, the calculation^ by 
the formula Eq. @ well explains the experimental re- 
sults quantitatively, giving some credit that the AHE is 
mostly of intrinsic origin rather than extrinsic origins, 
e.g., skew scattering and/or side-jump mechanism. How- 
ever in the presence of the time-reversal symmetry, the 
d.c. transverse conductivity a xy vanishes, and the topo- 
logical structure of the Bloch wavefunctions has not been 
systematically studied in the context of transport the- 
ory. As we will show below, an even more beautiful and 
nontrivial quantum topological structure is hidden in the 
valence-band structure in the paramagnetic state, which 
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is analogous to the fermionic quasi-particles in the SO (5) 
theory 2 ^. This is also motivated by the recent work by 
one of the present authors on the generalization of the 
quantum Hall effect into four dimensions in terms of the 
SO(5) symmetry 2 ^. In this paper, we shall show that 
the SO (5) group structure of the S — 3/2 Bloch states 
provides a natural description of the non-abelian SU (2) 
holonomy and its curvature in momentum space. This 
gauge structure underlies the dissipationless, topological 
spin current in hole-doped semiconductors. 

In the presence of the spin-orbit interaction, the con- 
ventionally defined total spin operator is not conserved, 
and it is nontrivial to define the spin current in this case. 
Our formalism resolves this issue by discovering con- 
served quantities in the Luttinger Hamiltonian (0J and by 
defining associated conserved spin currents. These quan- 
tities have clear physical and geometric meanings. The 
exact quantum calculation of the conserved spin Hall con- 
ductivity a s is performed in terms of the Kubo formula, 
and the results can be expressed entirely in terms of the 
non-abelian gauge curvature in momentum space. Our 
fully quantum mechanical results identify the quantum 
correction to the previous semiclassical result 2 in terms 
of the wave-packet formalism. 

The plan of this paper is as follows. In section [nj we 
reformulate the Luttinger Hamiltonian in terms of the 
SO (5) algebra and give the definition of the conserved 
spin current. Based on this definition, the calculation of 
the spin Hall conductivity in terms of Kubo formula is 
presented in section lTTTI where the geometrical meaning is 
stressed, and the comparison with the previous semiclas- 
sical result is given. Section Hvl is devoted to conclusions. 
Throughout the paper, we take the unit h = 1, e = 1. 



II. DEFINITION OF THE SPIN CURRENT 



From this projector form of the Hamiltonian, we see that 
the LH and the HH bands are each two-fold degenerate, 
and there is an SU(2) rotation symmetry acting on each 
band. Combining the LH and the HH bands, there is an 
SO(4)=SU(2)xSU(2) symmetry at every k point. In this 
section, we shall develop the mathematical framework in 
which this SO (4) symmetry is made manifest, and this 
symmetry is used to define the conserved spin current. 
Since P H and P L depend on k, the quantization axis for 
each SU(2) varies as a function of k. When k is adiabat- 
ically changed along a closed circuit, the fermionic wave 
function in general does not return to itself; in fact, the fi- 
nal wave function is related to the starting wave function 
by an SU(2) transformation within each band. Therefore, 
this problem is a natural generalization of Berry's U(l) 
phased to the case of SU(2) holonomy 22 i 25 . 26 i 27 i 28 i 29 . 30 . 
In particular, Demler and ZhangS developed a formal- 
ism of the SU(2) non-abelian holonomy in terms of the 
SO (5) Clifford algebra, which we shall adopt through- 
out this paper. Upon expanding the (kiS 1 ) 2 term in the 
Hamiltonian we obtain a product of two quadratic 
forms, one of the form Cakikj and another of the form 
£* J S' l iS y , where is a symmetric matrix. A 3x3 sym- 
metric matrix can be further decomposed into one trace 
and five traceless parts. The trace part has the same 
form as the first term in the Hamiltonian and can- 
cels the 72 contribution by construction. The remaining 
five traceless symmetric (£ l J — = 0) combination 

of S t S : > can be identified with the Clifford algebra of 
the Dirac T matrices, with the identification 



r a =e'{^,^'}, {T a ,T b } = 2S ab 



(11) 



The explicit forms of S -1 , T a and £ l J are given in Appendix 
lAl In terms of the F matrices, the Luttinger Hamiltonian 
(J3J takes the elegant form 



The Luttinger Hamiltonian JIJ has two eigenvalues: 

n /n 7l+ 2 72, 2 TP (U\ 71 - 2 72 2 , . 

E L{k) = — k , E H (k) = — k (7) 
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corresponding to the light-hole (LH) and the heavy-hole 
(HH) bands. Each eigenvalues are doubly degenerate, 
due to the Kramers theorem based on the time-reversal 
symmetry. For a fixed value of k, let P L (k) denote a 
projection onto the two-dimensional subspace of states 
of the LH band. We also define P H (k) similarly. These 
operators are written as 

9 
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2V k-S) 2 , P H = 1-P L (8) 



They obviously satisfy 

P L P H = = P H P L , {P L ) 2 =P L , (P H ) 2 =P H (9) 

In terms of these projectors, the Luttinger Hamiltonian 
can be expressed as 

H = J2(E H WP H (k)+E L (k)P L (k)). (10) 



H(k) = e(k) + ^d a T a , 
m 



(12) 



where 



<k) = ^k\ d a (k) 



di = 
di = 
d 5 = 



~3£,a ' kikj, 



-V3k v k z , d 2 = —V3k x k z , d 3 = -V3k x k v 
~~^~(k x ~ k y ), 



(13) 



We recognize that the vector components of d a (k) are 
nothing but the five d-wave combinations in the k space. 
The five-dimensional vector d provides a mapping from 
the three-dimensional k space to the five-dimensional 
d space (Fig. Since the Luttinger Hamiltonian de- 
pends on k only through d(k), we can perform all cal- 
culations in the 5D d space, and finally project back 
onto the 3D k space. This formalism enables a uni- 
fied treatment for the anisotropic Luttinger Hamiltonian, 
and more importantly, reveals the deep connection to the 
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four-dimensional quantum Hall effect (4DQHE)2^. Here 
and henceforth we adopt the convention that indices ap- 
pearing twice are summed over. 




FIG. 2: Mapping from the three-dimensional k vector space to 
the five- dimensional d vector space. The gauge structures of 
the spin current are associated with the magnetic monopoles 
located at the origins of the k and d spaces. 

The eigenvalues of (T2J are E L (k) = e(k) + ^d(k) and 
E H (k) = e(k)-£d(k), where d(k) = ^UT a = k 2 . They 
are of course the same as Q. In terms of the T matrices, 
the projection operators JHJ can be expressed as 

P L = \{l + d a T a ) , P H = l -(l-d a T a ) 1 (14) 

where d a = d a /d. The T matrices are convenient for sub- 
sequent calculations, since a product of any number of T 
matrices can be easily reduced to a linear combination of 
1, T a and r afc = j-[T a ,T b }. The five T a matrices contain 
the most general quadratic terms of the spin operator S l , 
while the ten T ab matrices contain both the three spin op- 
erators S l and the seven cubic, symmetric and traccless 
combinations of spin operators of the form S l S : >S k , as 
discussed in Appendix 1X1 These ten spin operators are 
generated under the Heisenberg time evolutions of the 
single spin operators, and it is natural to group them all 
into a unified object. For 72 = 0, r afc commutes with the 
Hamiltonian and generates an SO (5) symmetry group of 
the Hamiltonian. (In fact, the Hamiltonian has a higher, 
SU(4) symmetry in this case). For a given k, a fixed 
d vector singles out a particular direction in the five- 
dimensional d vector space, and the second term in <|12|) 
breaks the SO(5) symmetry to an SO(4) symmetry. This 
is nothing but SO (4)= SU(2)xSU(2) symmetry which we 
discussed earlier. In this way, we see that the SO (5) for- 
malism gives an elegant geometric interpretation of the 
SU(2)xSU(2) symmetry of the LH and the HH bands. 
It is in fact a subgroup of SO (5) rotation which keeps 
a fixed d vector invariant. As we shall see later in the 
paper, and in Appendix IU1 the 3D monopole structure in 
k space can be best understood in terms of the monopole 
structure in the 5D d space. 

To find the conserved quantities, let us define the con- 
served spin density explicitly as 

p ab (p, q) = cJ +?i/1 P o6)0d (p, q)r^C p _|, IO (15) 

where P a b,cd = P c d,ab = -Pba,cd = -Pab,dc- For the 
conservation of the spin, we require that — ip ab (p,q) = 
[H, p ab (p, q)] is proportional to q for small |q|. This is 



realized by imposing a condition on P a b,cd as 

Pab,cd(P> q) (dc ( P + |) + 4 (p - |) ) = (16) 

or equivalently, 

iV*(p>q)(d»(p + §) +d «(p-f)) =°- 

From these relations it follows that 

P ab (p, q) (d a (p + I ) + d a (p - I)) = 0. (18) 

There are five such linear equations, but only four of 
them are linearly independent because of the antisym- 
metry of p ab . Originally, p ab has 10 degrees of freedom, 
subtracting 4 constraints gives the remaining 6 degrees 
of freedom, exactly the same as the number of genera- 
tors in SO(4) algebra. Therefore, the projection opera- 
tor P a b,cd projects the full SO(5) symmetry generators 
into the SO(4) subspace which is orthogonal to a given 
direction of d. 

In the limit q = 0, l|17fl is satisfied by 

Pab,cd = £ (j>ae§bd — 5 a dSbc + 5 a ddbdc 

- Sbddadc - S ac dbdd + Sbcdadd^j ■ (19) 

It satisfies P a b,cd p cd,ef = P a b,ef, implying that it is prop- 
erly normalized as a projection. It is interesting to note 
that P a b,cd can also be expressed as 

Pab,cd — faba'b' fa'b'cdi ] abed — ~^^abcdede- (^0) 

Inserting l|19|) into the spin density l|15fl . we obtain 

p ab ^E^^kXr-Vck, 

k 

= E 4 ( ra " - \\ d * Tb - ^ r "' d f Y f ]) c k „(21) 

k \ / [iV 

Because [r ab , d e Y e ] = 2i(d a T b - d b T a ), we get 

p ab = E c L ( rab - \ti rab > d e nJ f rf]) c k „ 

k V / nv 

= E C L {P L T ab P L +p H r ab P H ) u c ku . (22) 

k 

Thus it corresponds to projecting out the inter-band ma- 
trix elements of F"' J . The conservation of p" b becomes 
manifest in H22J) . because the Hamiltonian is diagonal in 
each subspace, i.e. the LH or the HH band. 

The equation of continuity determines the uniform spin 
current to be 




pL T abpL + pH T abpH}^ ^ ^3) 
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To connect this spin current with the physical spin cur- 
rent in Ref. 0, we define a tensor rf a b by S l — rf ab T ab = 
5j7?ab[r a } r 6 ]. Explicit forms and properties of if ab are 
summarized in Appendix 1X1 By contracting with \rQ b , 
the conserved spin takes the form 



<?' — Ln l n ab 
3(c) — o^abP 



1 

3' 
3 V 



(P L S l P L + P H S l P H ) v c^ v . (24) 



The subscript ( c ) denotes the fact that this spin current 
is conserved. Here, we inserted a factor of ^, because 
in the LH and HH bands (i.e. S = 3/2 subspacc), the 
expectation value of the spin angular momentum is one- 
third of that of the total angular momentum S k . Thus, 
Eq. (|24|l corresponds to neglecting interband matrix ele- 
ments of the spin angular momentum. In a matrix form, 
the corresponding conserved spin current is 



ji _ 1 1 J ®H pL cl p l , pH ni p h\ 
J W~3 lldki' F SF +F bF ) 



(25) 



III. KUBO FORMULA CALCULATION OF THE 
SPIN CURRENT 



and G(k, iu) n ) is the Matsubara Green's function, given 

in pn) . 

In the clean case, the summation over w„ can be cal- 
culated by a contour integral. In the trace operation in 
the above equation, only the terms of products of four or 
five r matrices are nonzero, and the result is 



V 



d{uo 2 - 4-f%d 2 /m 2 ) 



\^-tijkk 2 kik k - e l j k k 2 k k ki j , (2!)) 



where ez,(k) = £x(k) — n and e#(k) = En(k) — H are 
one-particle energies for the two bands, measured form 
the chemical potential /i. 

Therefore, in the static limit the linear response is 
given by 

cr', = hm — J - 

J LJ^O -IUJ 



1 n F (e L ) - n F (e H ) ( 71 
3^4- fc 2 £yi V272 



-iei«(fcf - Ap) {— + 1 
6tt 2 3 V f Fy V 2 72 



(30) 



A. Difficulties with the conventional definition of 
the non-conserved spin current 

In order to calculate the spin current response based 
on Kubo formula, we should first define the "spin current 
operator" . The conventional definition of the spin cur- 
rent, with spin along the / axis flowing along the i axis, 
is given by 



J 



1 1 f dH 

"32 \dk~i 



(26) 



Because the spin S l is not conserved, J\ does not sat- 
isfy the equation of continuity without any source term. 
Before presenting the full calculation based on the con- 
served spin current discussed in the previous section, wc 
first calculate the linear response of this non-conserved 
spin current to the applied electric field and then com- 
ment on its difficulties. The Kubo formula gives 

1 r - 

Q\Svm)=-T7 (TJl(u)Jj)e* UmU du 



VJ 

= yoJ2 tx ( J i G ( k > *K + v m ))JjG(k, iuj n )) (27) 

k,n 

where v m — 2irm/[3 (m: integer), w n = (2n + l)ir/P 
(n: integer), (3 = l/fc^T, T in l|27() represents the time- 
ordering, 



j. - V c f [ — + V ^r h 



(28) 



This result a\j in Eq. I|3UI) does not vanish in the limit of 
72 — * 0, i.e., the absence of the spin-orbit coupling. It is 
not a contradiction, because the two limits 72 — > and 
oj — > cannot be exchanged in Eq. (|29|l . and Eq. (|30|l is 
the one which is valid in the d.c. limit, when —k 2 ^> u>. 
We have learned that Hu, Bcrnevig and Wu have also 
obtained a similar result independently^!. We note that 
this result l|3UI) is reproduced by wave-packet dynamics 
inrcf.Hi 

The conventional definition of the spin current (|2(j|| is 
physically admissible, as is usually adopted. However, 
its mathematical meaning as a "current" is ill-defined. 
A "current" is always associated with a corresponding 
conserved quantity. A "current" is then defined by using 
the Noether's theorem, or equivalently, by the equation of 
continuity. Since the conventionally defined spin current 
is not conserved for the Luttinger Hamiltonian due to the 
spin-orbit coupling, we shall use the the conserved spin 
current constructed in the previous chapter. 

There are also physical reasons to take this conserved 
spin current. Generally speaking, there must be some 
reason for a quantity to have slow dynamics and to con- 
tribute to the low frequency response. One is a conser- 
vation law and the other is a critical slowing down. In 
the present context, the latter is irrelevant and we need 
to look for a conserved current as we have done in the 
preceding chapter. When we separate the spin into the 
conserved and the nonconserved parts, the nonconserved 
part 



J i(n) 



1 



(P L S l P H + P H S l P L 



(31) 
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has an oscillating factor in time e ±l i E L-E H )t m 
Heisenberg picture. Its frequency is E^ — Eh and is nom- 
inally 0.1-1 eV or 1-10 fsec. As we are observing spins 
averaged over the time-scale much longer than l-10fscc, 
the only remaining part is the conserved part. Thus, in 
addition to mathematical soundness, the conserved part 
of the spin current automatically takes into account this 
averaging over time. In the next section we shall calcu- 
late the d.c. response in terms of the conserved part of 
the spin current. 



Kubo formula calculation for the conserved spin 
current 



In contrast with the previous approach, the approach 
using conserved spin Sl c \ gives well-defined and conserved 
spin current (|25|l . This approach is equivalent to neglect 
interband matrix elements of spin operators S l , as seen 
from l)24[l. This is justified in calculation of spin current 
because of the following reason. Let us consider the prob- 
lem in a semiclassical way. Two wave-packets in different 
bands are moving with different velocities, and they will 
move apart inside the sample. Meanwhile, in the sam- 
ple there are sources causing decoherence between wave- 
packets, e.g. inelastic scattering. This decoherence effect 
smears out the interband matrix elements. Therefore, in 
the measurement of the spin current, what is measured 
is only an intraband matrix element of spin carried by a 
hole coming out of the sample. Thus in the measurement 
of the spin current, we should consider only the intraband 
matrix element of S L . This is in contrast with calculation 
of susceptibility, where intraband matrix elements of spin 
gives significant contributions. 

By applying the electric field, this (conserved) spin cur- 
rent is induced by spin-orbit coupling. Let us calculate 
this linear response according to Kubo formula. Hence 
we shall calculate 



$(iv m ) = -- J {TJ?{u)J 3 y—du 

))J i G(k,iw„)).(32) 



V(i 



k.n 



By evaluating the summation over lo„ and taking the 
trace as presented in Appendix El we get 



(zi/„ 



16^ m /72\ 2 y^ n F (e L ) - n F {e H ) l2riab 
\m/ 4^ (i%) 2 — ^\d? jm 



V 



where 



(33) 



a6 _ 1 ddddde 



(34) 



is a purely geometric tensor. In the static limit we have, 

ab y 

cr 4 A = lim — 

= ±£(n L (k)-n H (k))G# 



V 



5> L (k) - n H (k))(F^ ab - F^), (35) 



where Ul = ^f(cl), n H = ^F(eff) ar e the Fermi func- 
tions of the LH and the HH bands. Here F^' ab and F?f ' ab 
are non-abelian gauge field strengths, i.e. curvature of 
the gauge field in the LH and the HH bands, and their 
definition and formulae are given in Appendix [O In con- 
trast to the result of the non-conserved spin current, the 
conductivity of the conserved spin current (|35f) is express- 
ible in terms of purely geometric quantities. Here we note 
that G° h as given in l|34l) is similar to the 9 term in the 
(l+l)-dimensional 0(3) nonlinear cr-model, which takes 
the form of 



drij drik 
dk a dkp 



(36) 



In fact, (|34|l describes the mapping of an area form 
from the three-dimensional (R 3 ) k space to the five- 
dimensional (i? 5 ) d n (k) space . An area element on R 3 
has 3 orientations dki A dkj , while an area element on R 5 
has 10 orientations, d(d a ) A d(db). Our formula describes 
the Jacobian of the area map. Out of the 10 possible ori- 
entations of an area form in R 5 , the f a bcd — \(-abcded e 
tensor in (|34|l selects 6 orientations which are locally 
transverse to d a . Geometric properties of the G° b tensor 
are further summarized in AppedixICl 

By substituting the formula l|C17|) for , we get 



(37) 



By contracting with ^rj l ab , the linear response of the cor- 
responding current is 



J = 



T) 1 n ah 



Gtt 



2 c y 



l{kp kf 



(38) 



where we used l|A33(l in Appendix El In contrast to the 
result (|30l) of the non-conserved spin current, the con- 
ductivity for the conserved spin current (|38() vanishes in 
the d.c. limit when the spin-orbit coupling 72 vanishes. 



C. Spectral representation of the response function 
in terms of the non-abelian gauge field 

The Kubo formula result for the conserved spin current 
obtained in the previous section can also be obtained by 
the spectral representation of the response function in 
terms of the eigenstates of the Hamiltonian. This treat- 
ment is similar to the one in quantum Hall effect by Thou- 
less et al^. By expressing the Kubo formula in terms of 
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the eigenstates, we can directly obtain the spin Hall con- 
ductivity in terms of the curvature F^ of the non-abelian 
gauge field for each band. 

Inserting a set of complete eigenstates into <|32|) . we 
obtain 



rv*(i ^- 1 V- ( {uLk\jf\pHk){pHk\J 3 \aLk) 



— — — — J - (n H - n L ), (39) 

-2 72 a/m + iv m J 



where |aLk) and \f3Hk) (a — 1,2, f3 — 1, 2) are the pe- 
riodic part of the Bloch wavefunction with wavenumvcr 
k in the LH and the HH bands, respectively. By substi- 
tuting 



j _ ®H jab _ 1 J ®H pL r ab pL , pH r ab pH 



(40) 



we get 



r" b = lim 



Therefore, in the d.c. limit 

QffM 

->0 — iuj 
/ d(fJLk) 
\ dk t 
, / d((3Lk) 



p 



H 



I 

2V 
d{aLk) 



dkj 



( /dffiHk) 
V \ dk 



P 



P l 



H 



d(aLk) \ 
dh / 
aLk|r ab |/3Lk) 
d(aHk) ^ 



/ d((3Hk) 
\ dk 



P L 



dk. 
d(aHk) 



dkj 



■(aHk\T ab \PHk)] 



(44) 



This formula can be expressed with the field strength Fij 
of the SU(2) gauge field for each band. We define the 
gauge field for the LH band as 



i ( aLk 



d 



PLkj , (45) 



SfclO = Wr £ (»* - 



2V 



n L ) 



a,/3,7,k 



fWLkm | 7 ffk)( 7 ffk||f|aLk) 



2 7 2<i/m + if,, 



(/?Lk|f£| 7 ffk)( 7 Hk|f£|aLk) 



2 7 2<i/m — if,, 



(aLk|r Qb |/?Lk) 



'C9Hk|^| 7 Lk)( 7 Lk|fg|aHk) 



\ 2 72 d/m + ii/,„ 

(/?gk|ff | 7 Lk)( 7 Lk|f||affk) 
2 72 (i/m— if m 
•(ai/k|r afc |/3i?k)] . 



(41) 



and similarly for A^. The corresponding field strength 
is 



J7 L — L 

13 dh 



dA[ 
dkj 



(46) 



and F.h , respectively. While in this definition Af is a 
2x2 matrix, it can be embedded into 4x4 matrix by 
identifying it with \aLk)(Af) a p(f3Lk\. We use the same 
notation Af to denote the 4x4 matrix defined in this 
way. The 4x4 matrices Af , F?j, and F^ are defined 
similarly. They can be expressed as linear combinations 
of T ab as 



zpL rpL.aby^ab zpH rpH,ab-pab 

i'j ~ ij ' U ~ ij ' 



(47) 



Then the resulting form of the spin Hall conductivity is 
obtained as 



It can be checked that Qfj(iv m = 0) = 0. Here we shall 
use the Feynman-Hellman theorem; because H\jHk) = 
Eu\jHk) implies 



it follows that 



d^Hk) \ _ dE H 
dk; / dh 



IjH^+Eh 



d(lHk) \ 

dh / 
(42) 



[PLk 



dH 



jHkJ 



dk, 

2 l2 d(k) / d(PLk) 
\ dh 



m 



2 72 d(k) 



7 ffk 



(3Lk 



d(lHk) \ 
dh J 



(43) 



ab 4 

•i =v 



J2(^-n L )(-Fh ab + F^ ah 



(48) 



in exact agreement with <|35|) . By contracting with \rf ah 
as in we get 



7 y(c) 



o r lab a ij 



^rfab J2( UH ~ UL "> (~ F > 



pL.a6 | ttiH, ab 
ij + t ij 



Fh 



(49) 



in exact agreement with (13811 . 



8 



D. Semiclassical limit 

The above result can be written as correlation func- 
tions in a real-time formalism; 

1 f°° 

4(c) = ^ I dt e l{uJ+U)t i{Mt), S\ c) }, J^e-P", 

(50) 

where Z = tre " is the partition function of the equi- 
librium. This quantity does not change if we replace S 1 ^ 

defined in 124|) by S'^ = Xk l , which follows from the fact 
that the helicity is a conserved quantum number. 

In a semiclassical (sc) approximation, one treats the 
spin S'As as a classical variable, commuting with the cur- 
rent Jj. Under this approximation, one obtains 

1 f°° 

(51) 

where we used the fact that S'^ commutes with H . Di- 
rect computation of this correlation function leads to the 
semiclassical result 

4(c) ( SC ) 

= ^7 E (n H txF?P H S l P H + n L trF£P L S l P L ) 

k 

= ^2^7(3^-^). (52) 

which agrees exactly with the semiclassical results^ based 
on the wave-packet equation of motion. The noncommu- 
tativity between the quantum spin and current operators 
contained in l|5UI) leads to a quantum correction 

A 4(c) = 4(c) - 4(c)( SC ) = -J^2 e ^( k F + *£)-(53) 

to the semiclassical result (|52|) . 

We would like to stress that this difference arises from 
the definition of spin current. In 149|) . we defined the spin 
current as an anticommutator between velocity and the 
spin as i|25f) . This definition of spin current amounts to 
taking the spin as a quantum average between the initial 
state and the intermediate state in the Kubo formula, as 
can be seen from Eq. 141|) . On the other hand, the semi- 
classical result 2 corresponds to taking the spin as that 
of the initial state. In this semiclassical formalism, the 
wave-packets with different helicities have the opposite 
transverse velocities with respect to the external electric 
field. 



IV. CONCLUSIONS AND DISCUSSIONS 

In the present paper, we studied the spin Hall effect in 
hole-doped semiconductors such as Ge and GaAs. The 
four valence bands, which are made out of p-orbitals with 



the spin-orbit interaction, consists of the doubly degen- 
erate heavy-hole band and light-hole band. (When we 
assume the inversion symmetry, the Kramers theorem re- 
quires at least double degeneracy at each k-point.) These 
two bands touch at the T-point. The effective Hamilto- 
nian describing these valence bands, so-called the Lut- 
tinger Hamiltonian, has a beautiful mathematical struc- 
ture described by the SO (5) Clifford algebra. At a given 
momentum k, the spin-orbit coupling singles out a fixed 
direction in the five-dimensional space of the d vectors, 
and breaks the symmetry down to SO(4)=SU(2)xSU(2). 
This symmetry property can be used to define conserved 
spin currents in both the LH and the HH bands. The 
quantum response of the conserved spin current can be 
calculated exactly within the Kubo formalism, and the 
result is summarized in Eq. 149|) . This result can be ex- 
pressed in terms of purely geometric quantities, or equiv- 
alently, in terms of the non-abelian Yang monopole field 
strength, defined in the five-dimensional space of the d 
vectors. This result also establishes the deep connec- 
tion between the spin current in the Luttinger model and 
the 4DQHE model of Zhang and Hu 2 ^, which also uses 
the Yang monopole as the non-abelian background gauge 
field. In the former case, the Yang monopole is defined in 
momentum space over the space of the five-dimensional 
d vectors, while in the latter case, the Yang monopole is 
defined in the real space. Magnetic monopole structure 
in the five dimensional momentum space has also been 
discussed by VolovikS. 

Our fully quantum mechanical results are compared 
with previous semiclassical one (1521 , and a quantum cor- 
rection due to the entanglement of spin and velocity is 
identified. The quantum correction can be traced to the 
non-commutativity and entanglement between the spin 
and the current operator. In physical systems where 
this entanglement is destoyed by some decoherence mech- 
anisms, the semiclassical result might be realized. In 
Ref. |3j, Culcer et al. developed a wavepacket formal- 
ism, and discussed the difference between our semiclas- 
sical result^ and the Kubo-formula result 1)31) [I using the 
conventional definition of the spin current. They incor- 
porated the nonzero correlation between spin and veloc- 
ity into a "spin dipole" and "torque moment" terms in 
their wavepacket formalism, and reproduced the Kubo- 
formula result Eq. I|30[> after also including a first-order 
field correction to the wavepacket spin. 

In the calculations of the spin current presented in 
this paper, we assumed an absence of impurities. On 
the other hand, we have also done a calculation includ- 
ing a scattering by randomly-distributed impurities. By 
assuming that the scattering potential is isotropic and 
accompanies no spin-flip, we calculated the spin current 
within the Born approximation and the ladder approxi- 
mation for the vertex correction. The self-energy obtains 
a finite imaginary part h/2r as usual, where r is a life- 
time. The vertex correction, on the other hand, vanishes 
due to the parity, namely because the Hamiltonian is an 
even function of k. Thus as far as the broadening of the 
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energy h/r is much smaller than the energy difference 
between two bands El — Eh, the spin current calculated 
in (|38|1 remains unchanged. The details of the calculation 
are involved and will be presented elsewhere. 

The dissipationless spin current discovered in recent 
theoretical works has many profound consequences both 
in fundamental science and in technological applications. 
However, in models investigated so far, there is still a fi- 
nite longitudinal charge conductivity and dissipation as- 
sociated with charge transport. A key objective along the 
current line of research is to identify spin-orbit coupled 
system with a gap in the electronic excitation spectrum, 
which might lead to quantized spin Hall effect, similar to 
the familiar quantized Hall effect. This exciting possibil- 
ity is suggested by the fact that is represented as the 
integral of the gauge curvature over the occupied states, 
and does not require the Fermi surface across which the 
particle- hole excitation occurs. 



APPENDIX A: V MATRICES AND RELATED 
IDENTITIES 



With the expressions for the S matrices 



I 2 



s z = 



s x = 



/ VI 

' 2 

V3 
2 



-IJ 



1 



2 

2 



\/3 / 
2 / 



S y 



we get 



V3, 



(Al) 



s 



x2 



V3 



1 - -<T Z ® <J Z 



5 

I' 



S* a = -^. a *®l-l a *®o* + ^ ) 
S z2 = <r z ®a z + ^, 



S x S y + S y S x = V3a y $ 1, 
S y S z + S z S y = V3a z <g> a y , 
S Z S X + S X S Z = V3a z ®a x , 



(A2) 

(A3) 
(A4) 



(A5) 
(A6) 
(A7) 



where a 1 (i — 1, 2, 3) are the Pauli matrices. Let us define 
the T matrices as 



r 1 
r 2 

r 5 



Since r Q r b + r f T Q = 25 ab , These five matrices generate 
the SO (5) Clifford algebra™. We shall define the traceless 
symmetric tensor by i.e. 



a z d 




= -^=(S y S z 

Vs 


+ S z S y ), 


(A8) 


a z d 


§ a x 


= -^=(S Z S 2 


+ S X S Z ), 


(A9) 


a y t 


51 = 


= -^(S X S V 


f S y S x ), 


(A10) 


a x C 


51 = 


= -^(s x2 - 


s y2 ), 


(All) 


a z d 


5 a z 


- b I' 




(A12) 



Explicitly they are written as 

1 ... 1 



(A13) 



2V3 



c xx _ (W _ 
?4 — — S4 — 



c xx _ syy _ 
S5 — S5 — 



J 2a/3 
1 



J ?3 



2vT 



2yT 

) S5 



G 



1 

3' 



and those obtained by = They form the vector 
representation of the SO (5) algebra, and are expressed 
as 4 x 4 Hermitian matices. When we define a repre- 
sentation in this space of 4 x 4 Hermitian matrices as 
T ab \A) = \[T ab ,A}), = A, T ab = i[r a ,r h ], it is shown 
to be a product of two four-dimensional spinor represen- 
tations of SO (5). This product of two spinor representa- 
tions can be classified into the irreducible representations 
of SO (5), and each irreducible representation is expressed 
as a product of the elements of the Clifford algebra. Thus 
4x4 = 1 + 5 + 10, where 4 is the spinor representation, 
1 is a trivial representation, 5 is a vector representa- 
tion spanned by T a , and 10 is an adjoint representation 
spanned by T ab . These matrices 1, T a , and r oh , span 
the space of 4x4 Hermitian matrices. Moreover, because 
pip2p3p4p5 _ _p a p roc j uc t f m ore than two T ma- 
trices can be written as a linear combination of 1, T a , 
and r afc . It is thus possible to write S l in terms of these 
matrices as 

S x = ^-1 <8> a x + -(a* ® a x + a y ® a y ) 

= ^r 15 -i(r 23 -r 14 ), (am) 

S y = ?—l <8> a y + ~{-a x <g> a y + a y ® a x ) 



.^_r 25 + I(r 13 + r 24 ), 



(A15) 



S* 2 = cr z ®l + il(g)a*=-r 34 -ir 12 . (A16) 
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These are used to calculate the correlation function in 
the Kubo formula. To formulate the problem in a co- 
variant fashion, we define rf ab as S % = rf ab T ab , where 
r ab = i[r a , T fc ] are generators of the SO (5) algebra, and 
rf ab = —t]l a . Nonzero components of rf ab are 



and 



x _ x 1 x _ 1 

Vl5 — ^23 — Vl4 — J) 



^25 



V3 j, _ 1 j, _ 1 



1 



1 



%4 — 2' ^ 12 — 4' 

and the ones obtained by 77* b — —r\ % ba - The ten T ab ma- 
trices contain both the three spin operators S l and seven 
cubic, symmetric and traceless combinations of the spin 
operators of the form S z S : *S k . These seven cubic opera- 
tors are 



x \3 



(s x ) 



7^3, 



.15 1 ' -pl4 



-r ll --r 



23 



(A17) 



( 5 y)3 = _IV^ r 25 + | r 24 + ^pl3 ; (Alg) 



/qz\3 _ 13 12 7 34 

{S X ,(S V ) 2 - {S z ) 2 } 



(A19) 

^WV 4 , (A20) 



(s z ) 2 - (s x ) 2 } = -^r 25 - ^r 24 , (A21) 



{s z ,(s x f - (s v ) 2 } = v^r 35 , 

gx gz 1 qz gx V j-i45 



(A22) 
(A23) 



There are several useful formulae for T a , which are 
used in the calculation in this paper: 

[r ab ,r c ] = 2i(5 ac r b -s bc r a ), (A24) 

{T ab ,T c } = e abcde T de , (A25) 

[r ab ,r crf ] = -2i(s bc v ad - s bd v ac 

-S ac T bd + S ad T bc ), (A26) 
{r afc , T cd } = 2e abcde T e + 25 ac S bd - 25 ad 8 bc ,(A27) 

tr(r a r") = 45 ab: (A28) 

tr(r a r fc r c ) = o, (A29) 

tr(rT<TT d ) = 4{5 ab 8 cd + S ad S bc - 6 ac S bd ),(A30) 

tr(r a r fc rT d r e ) = -u abcde . (A3i) 

By substituting S l — rf ab T ab into the commutation re- 
lation [S l , S-*] — ieijkS k , one can easily derive 



[if, rf 



(A32) 



where rf is a 5x5 matrix with components tf ab . In other 
words, the matrices —Airj 1 form the spin-2 representation 
of the SU(2) algebra. It can also be shown that 



vlbVL = -tr(?7y ) = -Si- 



(A33) 



(A34) 



Let us write down the formula for d a . We can easily 
check that 

C 3 C = ^(SikSji + SuS jk ) - -^M«- (A35) 
Then it follows that 

(gk i k j )r a = ggk i k j {s k , s 1 } 

= J(k-S) 2 -^fc 2 . (A36) 
Therefore, by substituting 



(k ■ S) 2 = -k 2 + 3Q' /■-, /.-,. I '" . (A37) 

into the Luttinger Hamiltonian @ and comparing it with 
(H2J, we get 



d a — 3^ k^kj 7 



(A38) 



in accordance with ljT3)l . This tensor ^ can be expressed 
in terms of r\ k d as calculated below. 



(k-S) 2 = 1 -{k t r lab T ab , kiffJT 1 } 



k 2 + kikj e abcde r] ab rf cct T e . 



By comparing with Eq. (|A37(I we get 

£,e — ~^ e abcde r lab T f cd - 

One can also check that 



(A39) 



(A40) 



(A41) 



APPENDIX B: DETAILS OF THE KUBO 
FORMULA CALCULATIONS 

The electron Green's function is written as 

_! ) 

iu n - H + n J ^ 
1 



G M „(k,iw„) 



(iu„ + fi — e(k)) 2 — 7|d 2 /m 2 



• (icj n + fx-e(k) + ^-d a (k)T a ) 

f(k, iu n ) ( 9 {k, tW„) + ^d a (k)T a ) . (Bl) 
V 777, / uv 
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In the clean limit, the Kubo formula calculation proceeds 
as follows 



1 



V J 



(TJi(u)Jj 



l du 



±- tr ( jfG(K i{u n + u m ))J,G(k, io; w )) 



yg^2f( k > + v m)) f ( k , iu n ) 



■tr 



0K% A TTl Qf\i% 



\g(k,i(u n + v m )) + — d g T 9 \ 
I m i 



de + l2ddh rh 



dh- 



oti 



72 



d,r l 



(B2) 



where we used (|A25f) . 

To evaluate the summation over uj n , we use a formula 

-p X! f( k > «K + v m ))f(k, iu) n )(Cg(k, iu n ) + D) 



to (=f^C + D) (n F (e L ) - n F {e H )) 
72 d(k) {{iv m f - 4 72 2 d(k) 2 /m 2 ) : 



(B3) 



where C, D are constants. By noting that the term pro- 
portional to <?(fc, iu! n + ii> m )g(k, iw n ) becomes zero in tak- 
ing the trace of the matrix, we have, 



2#(*fm) 



17 £ 



to n F (e L ) - n F (e H ) 



47|d 2 /m 2 ) 



■tr 



2 TO 



^1 4. Ti^.r' 1 

TO C?fcj 



2 



12 



d,r l 



(B4) 



The matrix inside the trace is a linear combination of 
products of two, three, four and five T matrices. By 
taking the trace, only the products of four and five T 
matrices survive. It is worth noting that the and -j^- 
gives no contribution; the former is because of P a b,cddd = 
and e c dghidgdi = 0, and the latter is due to d g T 9 diT l = 
d 2 . After some calculation it becomes, 



■16u„ 
V 



(H) v — 

Vto/ (iv. 



n F (e L ) - n F (e H ) 



(iv m ) 2 - <i%d 2 /m 



In the d.c. limit we have 



d 2 G?l 



(B5) 



(16(k ■ r,) 3 + fc 2 k ■ »7) o6 (n L - n H ){B6) 



where we substituted QC17|1 . Because of the spherical 
symmetry of the problem, the summation over k can be 
simplified further. By using identities 

k 6 k 

^2^(k)kikjk k ki 

k 

= ^(Sijdki + S lk 5 3l + 5 a S jk ) £ $(fc)fc 4 ,(B8) 

k 

where $(fc) is an arbitrary function of k = |k|, we can 
calculate as 

^^■ri)(n L -n H ) = lj:^L^ (B9) 



^^■/c;(k-?7) 3 (7i L -njj), 



L X] Z^VV + rfrfrf + r?rfrf){n L - n H ) 



15 ^ k 

k 



OAC\ 2^ 1.2 ^ 



240 ^ fc 

fe 



where we used (IA32I) (IA34I) . Hence 



(BIO) 



rr^-— lj c . V UL ~ UH - 4 n » e ,(k L -k H ^ 



k 



(Bll) 



APPENDIX C: MAGNETIC MONOPOLES IN 

d = 3 AND d = 5 

From Eq. (|12|l we see that the microscopic Hamiltonian 
depends on k only through the 5D vector d(k); therefore, 
it is natural to define the most general 5D gauge connec- 
tion in the d space, and then project the gauge connec- 
tion to the 3D k space. Let P L and P H the projections 
onto the LH and HH bands. These projections have the 
following properties; 

p L = \{i + d-v), p h = l(i-d-r) = i-p L , 



(p L y = p L , (p H y = p H , p"p l = o = p l p 

We can define the covariant gauge field strength, i.e. cur- 
vature F a b in terms of these projection operators as 



(CI) 



This gauge field is defined over the 5D d space, with 
spatial indices a, b = 1, 2, 3, 4, 5. It is a 4x 4 matrix, being 
a linear combination of the SO(5) Lie algebra matrices 





■QpL 


q p l- 




■QpH 


0P H } 


i 


dd a ' 


dd b _ 




dd a ' 


dd b 
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r . It can be explicitly evaluated as 



F n h — — 



1 



dd c v c dd d T d 

dd a ' dd b 



= ^(T ab + d c d b T ca -dJ a T eb ). (C2) 
It can also be written as 

Fab = Pab.cdX = fabef fefcd^ j (^^) 

where f a b g h is given in (|2"U1) . 

The gauge potential corresponding to the gauge field 
strength F a b is given by A a — — ^pdbT ab . This can be 
shown by explicit calculations, using the standard defini- 
tion 

F ^ = ^-^T +i ^ A ^ (C4) 
dd a od b 

From F a b we can define the dual field strength G a b by 



Gab — ^{F a b, rf c r C } = ^-^J fabcdT Cd , (C5) 

where we used Eqs. (|A25I) and i|C2l) . 

We now define the gauge field strength for each band 

as 



Fab 



fab 



-IP L 
-IP H 



-QpL 


QpL- 


dd a ' 


dd b 


QpH 


dP H 


dd a 


ddb 



It is easy to see that 

fab — -r a b + - t ab ) ^ah — f a b - fab- 



(C6) 
(C7) 

(C8) 



Since F a b and G a b are related to each other by a duality 
transformation 



fabcdG c d — F a bi fabcdFcd — G a bi 



(C9) 



F^ b and F^ are self-dual and anti-self-dual, in the sense 
that 



jabcdF c d — F a b J fabcdF l 



H 

cd 



^ab- 



(CIO) 



We can explicitly see that F^ b and F^ describes a gauge 
field strength with Yang monopole at d = 0. Let us 
define the two-form F L and F H as 

F L = \Ft b dda A dd b , F H = ^F%,dd a A dd b . (Cll) 

One can calculate that 

tr{F L A F L ) = -tr(F H A F H ) 
= -^dabcdeda ■ dd b A dd c A dd d A dd e . (C12) 



When this is integrated on a four-dimensional hypersur- 
face surrounding d = 0, it gives the second Chern num- 
ber multiplied by 87r 2 . Therefore F L and F H describe a 
gauge field with the Yang monopole at the origin, with 
its strength (i.e. the second Chern number) given by +1 
and —1, respectively^. 

Because of the projection operators P L and P H , F^ h 
and F^l can be expressed as SU(2) matrices operating 
within the LH and the HH bands respectively. In fact, we 
can see that they agree exactly with the conventional def- 
initions of the non-abelian holonomy or the SU(2) Berry 
connection. In the conventional definition, the SU(2) 
gauge field in the LH band as (A^) a = —i(aLk\ 3 ^ Lk ) 



dd a 

and its field strength is F^ b = d a A\ - d b A^+i[A^ A b ], 
where a, (3 — 1,2 characterize two eigenvectors forming 
the basis of the LH subspace. A^ and F^ can be defined 
in a similar way. The proof of the equivalence between 
the conventional definition and the definition <|C6I) can be 
seen in the following way, which is essentially the same 
as in Ref . El 



P 



■ dP L dP L 
dd a ddb 



L\2 



gpH QpL 

dd a dd b 



P- 



dP L ^ H 



dd a 



-P l 



dP L _ pL dP L {pH)2 dP L 



dd h 



dd a 



rtP L r)P H 
pL or p H or pL _ pL 



dd b 

dP L „ H dP L _ pL 



dd„ 

Ei aL > 



dd b 
d(aL 



dd a 



■P 1 



dd, 



dd a 



P 



H 



d((3L) 



dd b 



\ (f3L\, (C13) 



where \aL) = |aLk) and so forth. Then it follows that 
'dP L dP L 



dd a dd b 

= Ei« L > 

a, (3 

= Ei« L > 

a, S3 



— (a b) 



d{aL) 



dd a 
d{aL) 



8(fiL) 



dd b 



\ (pL\ - (a «- 6) 



dd a 



d{(3L) 



dd h 



(J3L\ 



8(aL) 



dd a 



- E \ aL ) 

a,/3,7 

— (a <-> b) 
i£|aL)(i^) tt/ ,</JL|, 

a, f3 



7 L ) ( 7 L 



(C14) 



which establishes the equivalence between (|C6ll and the 
conventional defini tion o f the gauge fields, for example, 
those used in Rcfs. 22 23. The equivalence between (|C7|) 
and the conventional definitions can be shown in a similar 
way. 

From these 5D monopole gauge fields, one can easily 
obtain the 3D monopole gauge fields by the pull-back 
mapping. For example, 



G 



dd a ddb_ 



« ~ dk, dk, Gab 



cdj^cd 



G%r 



(C15) 
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Substituting the definition of G a b as given in (|C5|) we see 
easily that G^ is given by Eq. (|3IJl . 

Calculation of G"? b and F°: b is straightforward but 
somewhat cumbersome. By using Mathematica, we ob- 
tain 



For example, we can take 

i + d 5 + ij: 4 a=1 r a5 d a 



u 



2(1 + 4) 



(C22) 



F$ = ie tfl *l(16(k • vf + 4^ 2 k • r}) ab , (C16) 

= ^e yi fci(16(k • r,) 3 + * 2 k • v)ab, (C17) 

where 77 — (ri x , r] y , rj y ) and ry 1 is regarded as a 5x5 spin- 
2 representation of the SU(2) Lie algebra, satisfying the 
commutation relation 1|A32|) . In these formulae, F a b and 
G a b are written in terms of T matrices. Alternatively, we 
can write them in terms of the spin matrices 5 1 ; 



(C18) 
(C19) 



1 L>A ~ I ) ( '.H-j..:- 



where A = k ■ S is the helicity matris^. Eq. (|C18|) has 
been obtained in Ref. I27t one can sh ow E q. fC19)l m 
the similar way. Equivalence between (|C18(I . (|C19|I and 
(|UT6|l . (|C17|l ca n be sh own b y sub stitut ing S l = ri ab T ab 
and using l|A27l) . From l|C18|) and (|CJ19|) . we get 



F H 



A 



A2 
2 



1 



c l _ T?L 



A 



3A^ 
2 



27 



£y7 



fc 3 ' 



(C20) 

As is expected, F" = for the LH band (A = ±1/2), and 
Ft = for the HH band (A = ±3/2), This is the field 
strength of the U(l) (Dirac) monopole with monopole 
strength ±3 for A = ±3/2 (HH band) and =f3 for A = 
±1/2 (LH band). 

Finally we would like to establish the exact equivalence 
between the gauge fields introduced above and the Yang- 
Mills instanton in Euclidean four-space^ or the Yang 
monopole gauge fields over the four-sphered. The proof 
essentially follows that of Jackiw and Rebbi 3 ^. The 2- 
form SO (5) gauge field on R 5 can be converted to SO (4) 
2-form gauge field on R 5 by gauge transformation U such 
that: 



By this gauge transformation, the gauge field A a and the 
field strength F a b are transformed to 

817 

U^A a U-iU^— = A a: (C23) 

Od a 

UjF ab U = F ab . (C24) 
These quantities A a and F a b are linear combinations of 
T mn (m,n = 1,2,3,4), belonging to the SO(4) algebra. 
Explicitly they are written as 

1 4 

Aa = -^ frE^ 1 " 6 (a =1,2,3,4), i 5 =0, 

2d(l + d 5 ) fri 



P ^ = -^T, dbTab (« = 1,2,3,4), 

b=l 



Fab = 



^{ Tab -ThP c{dbrac ~ daTbC) ) 



(0,6=1,2,3,4), 



u^d a r a u = r 5 



(C21) 



which are exactly the SO(4)=SU(2)xSU(2) gauge fields 
used in the context of 4DQHE2 3 .. 
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